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ABSTRACT 

Several families of irregular moons orbit the giant planets. These moons are thought to 
have been captured into planetocentric orbits after straying into a region in which the 
planet's gravitation dominates solar perturbations (the Hill sphere). This mechanism 
requires a source of dissipation, such as gas-drag, in order to make capture permanent. 
However, capture by gas-drag requires that particles remain inside the Hill sphere long 
enough for dissipation to be effective. Recently we have proposed that in the circu- 
lar restricted three-body problem particles may become caught up in 'sticky' chaotic 
layers which tends to prolong their sojourn within the planet's Hill sphere thereby as- 
sisting capture. Here we show that this mechanism survives perturbations due to the 
ellipticity of the planet's orbit. However, Monte Carlo simulations indicate that the 
planet's ability to capture moons decreases with increasing orbital eccentricity. At the 
actual Jupiter's orbital eccentricity, this effects in approximately an order of magni- 
tude lower capture probability than estimated in the circular model. Eccentricities of 
planetary orbits in the Solar System are moderate but this is not necessarily the case 
for extrasolar planets which typically have rather eccentric orbits. Therefore, our find- 
ings suggest that these extrasolar planets are unlikely to have substantial populations 
of irregular moons. 
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1 INTRODUCTION 

While the large regular moons of the giant planets follow 
almost circular, low inclination, prograde orbits, the aptly 
named irregular moons tend to do the opposite: i.e., they 
often have highly eccentric, high inclination orbits which 
may be retrograde or prograde. These moons, therefore, pro- 
vide actual observational examples of the complexities of 
three-body dynamics. Irregular moons are thought to have 
been captured during the early stages of the Solar System 
but the detailed mechanism has not been fully elucidated. 
Based on a study of the circular restricted three-body prob- 
lem we have recently proposed that chao s played a key role 
in the initial stages of capture fiAstakhov et al. 2003) ; in this 
mechanism - called chaos assisted capture (CAC) - parti- 
cles may initially become entangled in chaotic layers which 
separate directly scattering from regular (bound) regions of 
phase space. This temporary trapping serves to extend the 
lifetimes of particles within the Hill sphere thereby provid- 
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ing the breathing space necessary for relatively weak dissi- 
pative forces (e.g., gas-drag) to effect permanent capture. 
While this basic scenario may be responsible for the ini- 
tial formation of families of irregular moons, other factors 
may affect their subsequent long term evolution including 
collisions; gravitational perturbations from other planets; 
and deviations from circularity of the planet's orbit. Here 
we investigate the effect of orbital ellipticity of the planet 
on the capture mechanism using Monte Carlo simulations 
and phase space visualisations based on the Fast Lyapunov 
Indicator. Two representative star-planet binaries are con- 
sidered: the actual Sun- Jupiter system and a hypothetical 
extrasolar cousin of higher orbital eccentricity. We find that 
the chaos-assisted capture mechanism is robust to moderate 
ellipticities. 

Recent dis coveries of numerous i rregular moons at 
giant plan ets fGladman et al.l 1200 ll : ISheppard fc Jewitd 
120031 : iKavc laars et al. 20Qj|) have spurred the devel- 
opment of theories of their orig in, which is gener- 
ally thought to involve capture ([C arrubae^lJ j^OO^: 
l Yokovama et al J 120031: lAstakhov et all 120031: ICuk fc Burns! 
120041: iNesvornv. Beauee fc DonesI 120041: iNeto et alJ l2004f 
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iKavelaars et al.ll2004f) . These moons are considered to be 
primordial pieces of the Solar System captured in the am- 
ber of time. While there is, as yet, no consensus on the 
detailed capture mechanism of these minor bodies a number 
of recent studies have, nevertheless, attemped to simulate 
their post-capture evolution in order to explain salient as- 
pects of their orbits. For example, several attempts have 
recently been made to explain the observed clustering of ir- 
regular satellites as being the result either of the breakup 
of a larger parent bo dy or from catastrophi c collisions be- 
tween planetesimals I Nesvorny et al.| 12003: Cuk fc Burns! 



120041: iNesvornv et al.1 12OO4I : IKavelaars et al.l l2004^ . These 
hypot heses are suppor ted empirically by photometric sur- 
veys jGrav fc Holmanlt2004.'l which indicate that members 
of the clusters have similar surface colours suggesting that 
they may have had common progenitors. However, what- 
ever post-capture orbital evolution occurs within the Hill 
sphere there has to have been an initial capture mecha- 
nism efflcient enough to populate the Hill sphere with moon- 
lets in the first place. We have recen tly proposed one such 
mechanism - chaos-assisted capture iAstakhov et alJl2003fl 
- in which long-term, but temporary, trapping occurs in 
the 'sticky' ch aotic layers lying close to Kolmog orov-Arnold- 
Moser (KAM, iLichtenberg fc Liebermanlll992ll tori embed- 
ded within the planet's Hill sphere. If a particle gets entan- 
gled in one of these layers then even a moderate l evel of dissi- 
patio n (or, alternatively, slow planetary growth. iNeto et al] 
120041) may be sufficient to make capture permanent. It is 
likely that Comet-Shoemaker-Levy 9 was a recent example 
of an object trappe d in such a chaotic zone fo r the best part 
of the last century llChodas fc Yeoman JIQqH) . However, the 
absence of a large gas cloud at Jupiter makes contemporary 
permanent capture of such objects by gas-drag unfeasible. 



Qualitatively, in the CAC model the capture proba- 
bility of a planet is largely determined by the volume of 
phase space (at a given energy) taken up by chaotic zones 
within the Hill sphere. Simulations in the circular restricted 
three-body problem (CRTBP) show that capture becomes 
less probable at higher energies because most of the phase 
space is directly scattering and relatively few KAM tori 
survive; those that do are surrounded by relatively thin 
chaotic layers which ar e less effective at trapping intruders 
iAstakhov et alJ l2003l) . Factors that are not included 
in the CRTBP may also afi'ect the capture mechanism, 
e.g., the eccentricity of the parent planet's orbit. A more 
realistic extension of CAC beyond the CRTBP is, therefore, 
needed. As a first step we extend the mo del to the e lliptic 
restricted three-body problem (ERTBP. ISzebehelvl Il967t 
i Llibre fc Pifioll 119 90*: 'Benest' '200^; 'Pilat-Lohinger et al' 
I2OO3I: IVillac fc Sch ccrcs 2004 ; Mako fc Szcnkovits 2004) 
which allows us to model perturbations induced by 
deviations of the planet's orbit from circularity. While 
orbital eccentricities are relatively small, although not 
negligible, for the giant planets of the Solar System they 
may be appreciable fo r extrasolar planets (Schneider 1999.: 
Marcv fc Bulled l200d: IChiang. Fischer fc ThommesTfeoOs ' 



Mgrz^^^^^gi^gngghiUini 

Beaxige fc Michtchenkol |2 
2OOJ. For example. 



I2nn2t IColdreicli fc Sa.ril 1^ 
iMichtchenko fc Malhotral 

compare Jupiter's mean orbital 



eccentricity ^ e = 0.04839 to typical extrasolar planet 
eccentricities which lie in the range e 0.2 — 0.6 

iTremaine fc Zakamskalf2003t) . Up to now there have been 
more than 100 extrasolar planetary systems detected 
^ and an intriguing question is whether these massive 
planets might harbour moons jBarn'es' fc O'BrienI l2002l: 
iBurns fc CukI 120021: Iwmiamsll2003ll wh ich might be habit- 
able d WilUams. Kasting fc Wadd Il99i1) . As an exemplary 
extrasolar captor, we study a system similar in mass ratio 
to the idealized Sun- Jupiter binary but which follows a 
highly elliptic orbit. This provides a comparative estimate 
for the probability of capture and also suggests the ranges 
of energy and orbital inclination over which extrasolar 
irregular moons might be expected to exist. 

The paper is organized as follows: Section 2 intro- 
duces the Hamiltonian for the ERTBP which, in the limit 
of zero ellipticity, reduces to the CRTBP. We work in a 
coordinate system whose origin coincides with the planet 
and, in actual integrations, regularise the dynamics to 
deal with two-body collisions (|Aarsethll2003ll . Because the 
ERTBP is explicitly time dependent it is not possible to 
compute conventional Poincare surfaces of section (SOS, 
ILichtenberg fc LiebermanI Il992h even in the planar limit 
so as to visualise the structure of phase space. There- 
fore we use the notion of a Fa st Lyapunov Indicator (FLI, 
iFroeschle. Guzzo fc Leg'a ''2000'l to visualise the structure of 
phase space. Section 3 briefly discusses the CAC mechanism 
as applied to the CRTBP in order to facilitate comparison 
with the ERTBP. In particular, FLI on the surfaces of sec- 
tion are computed which can be compared directly with SOS 
in the CRTBP. This is done in Sec. 4 where Monte Ca rlo sim- 
ulatio ns of capture are performed. Unlike in Astakhov et alJ 
i2003h in these simulations we do not include dissipation 
and, instead, focus on the distributions of particles that can 
be temporarily trapped in chaotic zones for very long time 
periods. This avoids complications ass ociated with the bes t 
choice of dissipative force (see, e.g., ICxik fc Burn j l2'004^ . 
Finally, conclusions are contained in Sec. 5. 



2 HAMILTONIAN AND METHODS 

The CRTBP describes the dynamics of a test particle hav- 
ing infinitesimal mass and moving in the gravitational field 
of two massive bodies (the 'primaries' - e.g., a planet and a 
star) which revolve around their center of mass on a circular 
orbit. The equations of motion are, therefore, most naturally 
presented in a non-inertial coordinate system that rotates 
with the mean motion of the primaries ( Murrav & Dcrmot^ 
^9^. In the rotating coordinate system the positions of 
the primaries are fixed. When the planet's orbit is elliptic 
rather than circular a nonuniformly rotating-pulsating co- 
ordinate system is commonly used. These new coordinates 
have the felicitous property that, again, the positions of the 
primaries are fi xed; however the Hamiltonian is explicitly 
time-dependent JSzebehelvll 19671) . 



^ 'http : //nssdc . gsfc .nasa. gov/pla netary/planetf act .html | 
^ http://www.obspni.fr/plaiiets 
http : //exoplanets . org 
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2.1 Hamiltonian 

For our purposes it is most convenient to locate the origin at 
the planet (Fig. because angular momentum w ill be mea- 
sured with respect to the planet. Then, following ISzebehelvl 
il967f) andlEjibrc & Pinol ( 1990) , we obtain the planetocen- 
tric ERTBP Hamiltonian He after introducing an isotropi- 
cally pulsating length scale 



He ~ Ee 



;{{P'.c +yf + {Pv ~xf +p/ ■ 



+ ■ 



-(1- 



/(1 + ecos 



The semimajor axis of the orbit of the primaries Up has 
been scaled to unity and e is the eccentricity of the planet's 
orbit; p = mi /{mi + m^), where /i,mi and m2 are the 
reduced mass and masses of the planet and star, respectively 
[n = 9.5359 X 10"* for Sun- Jupiter). The true anomaly /, i.e. 
the planet's angular position measured fr om the pericente r, 
is related to the physical time t through l)Szebehelvlll967l ) 



dl _ (1 + e cos ff 
dt 



2)3/2 



(l-e^) 

The ( planet-centered) 
jAstakhov et al. 2003) 



(2) 



CRTBP Hamiltonian 



He^Ee=-{px +Py + Pz ) - {x Py ~ y Px) - 

1 - M M 



a/(1 -I- a;)2 + y'^ + z'^ y/x^ + y'^ + z^ 



(1 - n)x - -(1 - n), 



(3) 



is recovered when e = 0. In both cases (0 and @, 
hz = xpy — ypx is the z-component of angular momen- 
tum h — {hx, hy, hz) with respect to the planet. Note that 
the orbital inclination i — arccos /iz/|h| is invariant under 
isotropic pulsating rescaling from CRTBP to ERTBP. The 
orbit is said prograde if z < tt/2 {hz > 0) and retrograde 
otherwise. 

The Hamiltonian of the elliptic problem (0 is a periodic 
function of the true anomaly / (which plays the role of time) 
and, hence, generates a non-autonomous dynamical system. 
Unlike the circular problem, the ERTBP does not possess 
an energy integral and, evidently, no such useful guiding 
conce pt as a static zero- velocity surface jMurrav fc Dermotd 
I1999I) can be introduced. Furthermore, due to the extra di- 
mension associated with the explicit time dependence, con- 
struction of the SOS and, indeed, any visual analysis of 
phase space seems impossible even in the planar limit. 



2.2 Numerical methods 

We have performed numerical integrations in which fluxes of 
particles are simulated as having passed from heliocentric or- 
bits into the region surrounding the planet as defined by the 
Hill sphere. This region roughly corresponds to the region 
between the points labelled Li and L2 in Fig.0 In practice 
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Figure 1. Level curves of the zero velocity surface at high (a) 
and low (b) energies in the planetocentric planar CRTBP along 
with the five Lagrange equillibrium points (labelled Li - L5). The 
circle shaded grey is the Hill sphere centered on the planet (P). 
The star (S) is at (-1, 0). 



it is convenient to go to extended phase space by introduc- 
ing an additional pair of conjugate variab les - 'coordinate' 
/ an d 'momentum' pf = —H iLichtenberg fc LiebermanI 
Il992h for which the new Hamiltonian becomes conservative. 
In our numerical simulations we have used this method in 
combination with Levi-Civita (in 2D) a nd Kustaanheimo- 
Stiefel (in 3D) regularising techniques fetiefel fc Scheifel3 
119711 : lAarseth 2003) to avoid problems associated with two- 
body collisions jNagleilEooil . Numerical integ rations were 
done using a Bulirsch-Stoer adaptive integrator jPress et alJ 
tl999). 

One approach to visualising phase space structures in 
systems with greater than 2 degrees-of-freedom relies on 
computations of short-time Lyapunov exponents (LE) over 
sets of initial conditions of interest (for recent dcvcloi: 
merits and relevant applications sec. 
12000); S andor et alJ ll2 00n 
FCincotta. Giordano fc Simd 
and references therein, althou p^h the idea of using finite-time 
LE itself dates back at least to'Lor enzll965ll . As was demon- 
strated by Froeschlc et al. (2()00) various dynamical regimes 
(including resonances) can be distinguished by monitoring 
time profiles of a quantity called the Fast Lyapunov Indi- 
cator. Given an n-dimensional continuous-time dynamical 
system. 
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Figure 2. Poincare surfaces of section in the planar CRTPB at four energies (Jacobi constants Cj = —2Ec). Initial conditions were 
chosen randomly inside the Hill radius (here Rh is scaled to 1) and integrated with the cut-off time Tcut = 36000 years or until trajectories 
escaped. The SOS is the x — y plane with px = 0,y > 0. Points on the surface are coloured according to the sign of angular momentum 
hz as the trajectories penetrate the x — y plane (black, prograde with > 0; grey, retrograde with < 0). 



dx/dt = F(x, i),x = {xi,X2, ■■■,x„), (4) 

the F ast Lyapunov Indicator is defined as jFroeschle et alJ 
I2OOOI) 

FL/(x(0),v(0),t) = ln|v(t)|, (5) 

where v(t) is a solution of the system of variational equa- 
tions jTancredi. Sanchez fc R'oi3l200lfl 



dF 
9x 



(6) 



The complementary system (|SJ contains spatial deriva- 
tives which only aggravate the singularities in the restricted 
three body problem (RTBP) and so the use of regularisation 
becomes even more important. For this reason, all calcula- 
tions of FLI reported herein were made using regularised 
versions of @ and 0. 

As distinct from the largest Lyapuno v characteristic ex- 
ponent iLichtenberg fc Liebermanlll992l) 

v(t)| 



A — lim - In , , , , 
t-,00 t |v(0)| 



(7) 



which just tends to zero for any regular orbit, the more sen- 



sitive FLI can help discriminate between resonant and non- 
resonant regular orbits (Frocschle_et al. 2000). Although, in 
practice, detection of resonances may be tricky, since real 
differences in FLI for quasiperiodic and periodic orbits of 
the same island are not always huge. 



3 CHAOS ASSISTED CAPTURE 

In this section we briefly summarize the CAC mecha- 
nism described in I Astakhov et all (j200^ as applied to the 
CRTBP. Even though the basic mechanism applies in three- 
dimensions, it is simpler here to outline the general scheme 
in terms of the planar version of the C RTBP. Fig. [J shows 
the r elevant zero velocity surface (ZVS. lMurrav fc Dermoid 
ll_999|) together with the five Lagrange equilibrium points. 

Level curves of the ZVS, similar to a potential energy 
surface (PES), serve to limit the motion in the rotating 
frame and so define an energetically accessible region that 
may intersect the Hill sphere. However, unlike a PES, be- 
cause the ZVS is defined in a rotating frame, it is possible 
for energy maxima to be stable as is, in fact, the case for L4 
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Figure 3. A paradigmatic escaping 'sticky' trajectory (a) in the 
planar Sun-Jupiter CRTBP; its multiple returns (b) to the Hill 
sphere (shaded grey, Rh — 0.068) after heliocentric excursions; 
its radial distance from the planet r = ^ -\- xp- (c) and time 
profile of the Fast Lyapunov Indicator (d). The dotted vertical 
line indicates the first crossing of the Hill sphere (escape). 



and L5 a t which points Jupiter's Trojan asteroids are situ- 
ated Imirray & Dcrmott 1999). The Hill sphere (see Fig.Q 
roughly occupies the region between the saddle points L\ 
and Z/2 which we term the 'capture zone' with its radius 
being given by Rh = 0^/(^/3)^^^ where c tp is the planet's 
semimajor axis jMurrav fc Dermotd Il99^ . In the case of 
ERTBP the ZVS pulsates, which defines periodically time- 
dependent capture regions llMako fc SzenkovitdEool) . The 
two Lagrange saddle points L\ and L2 act as gateways be- 
t ween the Hill sph e re and heliocentric orbits. A key finding 
of lAstakhov et all (l2003h is that at energies close to (but 



above) the Lagrange points only prograde orbits can enter 
(or exit) the capture zone. At higher capture energies the 
distribution shifts to include both senses of /iz until, finally, 
retrograde capture becomes more likely. The statistics of in- 
clination distributions will, therefore, be expected to depend 
strongly on energy, i.e., how the curves of zero- velocity in- 
tersect the Hill sphere. Fig.|5|portrays the structure of phase 
space in the planar limit (z = Pz = 0) in a series of Poincare 
surfaces of section at four energies. At the lowest energy 
shown in Fig. |5Ja) many of the prograde orbits are chaotic 
whereas all the retrograde orbits are regular. Because incom- 
ing orbits cannot penetrate the regular KAM tori, prograde 
orbits must remain prograde while retrograde orbits cannot 
be captured nor can already bound retrograde orbits escape. 
After L2 has opened in Fig.|2Ib) the chaotic 'sea' of prograde 
orbits quickly 'evaporates' except for a 'sticky' layer of chaos 
which clings to the KAM tori. With increasing energy this 
front evolves from prograde to retrograde motion. Chaotic 
orbits close to the remaining tori can become trapped in 
almost regular orbits for very long times. In the presence 
of dissipation these chaotic orbits can be smoothly switched 
into the nearby KAM region and almost always preserve the 
sign of angular momentum. Thus, at low (high) energy per- 
manent capture is almost always into prograde (retrograde) 
orbits. 



4 RESULTS AND DISCUSSION 

In this section we describe the results of our simulations 
in the planar ERTBP using the FLl and our Monte Carlo 
simulations in the spatial ERTBP. 



4.1 Planar ERTBP and FLI 

We are primarily interested here in obtaining a qualitative 
picture of the volume and structure of phase space occu- 
pied by the chaotic layer in the ERTBP as compared to the 
CRTBP. For this we find a good diagnostic to be FLI{t) 
whose increase for chaotic orbits can usually be detected 
before, or no later than, a test particle finally escapes from 
the capture zone (Fig.|2Jl. These measurements, when made 
over the capture region (which contains permanently bound 
regular trajectories even at energies well above L\ and L2), 
provide an estimate of the number of orbits that could be 
captured as illustrated in Fig. 01 The number of chaotic or- 
bits inside the capture zone, detected by computing FLl, 
decreases with increasing ellipticity for moderate eccentrici- 
ties (Fig.^Ia)) signifying that the chaotic layers get weaker 
and, therefore, capture is expected to become less probable 
than it is at e = 0. 

Ideally the problem of separating the fractions of tem- 
porarily trapped and almost immediately escaping trajec- 
tories could be better quantified by partitioning the phase 
space into disjoint (e.g. 'inner'. Hill sphere, and 'outer', he- 
liocentric space) regions and computing the fluxes across the 
barriers between them. But, given the current state of phase 
space transport theory, this has not yet been shown possible 
in practice for essentially 3D problems such as the spatial 
CRTBP and ERTBP Hamiltonia ns. Interestingly enough, 
despite strong theoretical grounds iWiggins. Haller fc Mezig 
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Figure 4. Fraction of chaotic orbits within the Hill sphere 
that have FLI > 8 at the cut-ofT time Tcut = 200 years 
in the planar Sun- Jupiter ERTBP as a function of eccentric- 
ity. The phase space was sampled randomly with initial ener- 
gies, true anomalies and eccentricities taken also at random. For 
reference, shown in parentheses are the mean eccentricities of 
Jupiter, Saturn, Neptune, Uranus, Pluto, some extrasolar plan- 
etary systems (HD 209452, HD 169830, HD 73526, HD 108147, 
HD 74156b), and binary trans-Neptunian objects on elliptic orbits 
(2001 QC297,1998 WW31, 2001 QW322, 2001 QC298, 1999 TC36, 
1998 5Mi65). 



Il994f) . and exhausting attempts, efforts to describe quanti- 
tatively spatial three-body dynamics by constructing global 
invariant manifolds (iBelbruno 2004) . that would presumably 
contain all possible incoming and outgoing chaotic trajecto- 
ries through Li and L2 saddle points, have been not quite 
successful in approaching the 3D problem so far. Part of the 
reason is that multiple escapes and recurrences of high en- 
ergy trajectories to the Hill sphere (see example on Fig.|H^,b) 
make local manifolds near multidimensional saddles ex- 
tremely difficult to use as rigorous surfaces of no return. This 
is even more pro nounced for chaotic ionization of atomic Ry- 
dberg electrons llBrunello. Uzer fc Farrellvlll997l : iLee et all 
|200[|) which is closely related to RTBP dynamics. Recent 
progress in pursuit of manifolds for the spatial three-body 
problem is reported by 'Gomez et al." '2003*: l y iliac fc ScheeresI 
L2OO4; Waalkc ns. Burbanks fc Wiggin s 20ol 

To analyse the structure of phase space in ERTBP, 
we first computed FLI in the planar (2D) circular case 
{e = 0, z — pz =0), where direct com parison with surfaces 
of section (in the Hill limit, m < 1, ISimo fc Stuchill2000l: 



lAstakhov et all2003l) can be made. Figs.|^a-cl confirm that 
phase space visualisation via computing short-time FLI 
works well in the planar CRTBP, reproducing correctly all 
the relevant features visible in the SOS shown in Fig. |21 In 
particular, the chaotic layer, and its evolution with increas- 
ing energy, can be easily identified by high values of FLI 
(shown in yellow). Note also that FLI measurements make 
sense even for relatively short time integrations, and are 
more economical than constructing the corresponding SOS 
(compare cut-off times given in the Fig. |3 and Fig. |5| cap- 
tions). As for the sensitivity, it is not clear how reliable the 
numerical distinctions (red-purple-blue) between the FLI for 
quasiperiodic and resonant orbits are, but this is irrelevant 
for our current purposes. 

In the planar limit, initial conditions were generated 
randomly within the Hill radius on the surface of section 
(see Fig. (^caption) and assuming that, initially, / = n/2 in 
which case the ERTBP initial conditions reduce to those of 
CRTBP. This guaranties that aU ERTBP initial conditions 
are generated with identical initial energies p/(0) — — i5e, 
and these are, in fact, true CRTBP energies Ec- The set- 
ting, thereby, allows for a direct comparison between the 
SOS of Fig. |5| and results obtained using the FLI. Due to 
the dimensionality of the ERTBP, relaxation of the above 
constraint, e.g. choosing the initial true anomalies at ran- 
dom, produced FLI maps with no distinct structure. On the 
contrary, consistency in initial conditions at fixed / = 7r/2 
helps track the smooth changes from CRTBP to ERTBP. 

As the planet's orbital eccentricity is increased to reach 
its actual value for Jupiter, FLI maps reveal a reduction 
of the density of orbits within the chaotic layers (yellow 
paterns) visible on Fig. Eld-f) at each of the correspond- 
ing energies shown. However, at this moderate eccentric- 
ity, the phase space structures responsible for CAC ('sticky' 
KAM tori surrounded by the chaotic layers) generally sur- 
vive any deviations introduced by ellipticity. This suggests 
robustness of the CAC mechanism with respect to actual 
ellipticities of the giant planets' orbits which lie well be- 
low e ~ 0.4 (see Fig. 2J- We further confirm this by direct 
Monte Carlo simulations of capture probability (Sec. 4.2) 
which reveal that CAC indeed survives weak ellipticity in 
the RTBP. An important implication for future studies can 
be drawn by noticing that the heliocentric orbits of Pluto 
(with Charon) as well as of some recently discovered bi- 
nary trans-Neptunian objects ^ are quite eccentric (see ex- 
amples on Fig. |1J, but not to the extent characteristic of 
extrasolars (although the relative orbit of the binary part- 
ners as distinct from the heliocentric orbit of the combine 
can be very eccentric). These binary objects can be viewed 
in the framework of the Hill approximation too, so the 
weakly elliptic chaos-assisted mutual capture, probably sta- 
bilized by the fast exchange of energy and angular momen- 
tum with a fourth body, may have been involved in early 
stages of their formation ( Goldrcich. L rthwick fc Saril2002 : 
lAstakhov et aL,.2003: Funato et al. 200^. 

The picture at much higher eccentricities evolves to- 
wards significant distorsions of the phase space stuctures 
as visualised by FLI. In Fig. Elg-i) the KAM tori (red- 
blue islands) shrink, giving way to regions of scattering and 
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Figure 5. Colour coded Fast Lyapunov Indicator on the grid of isoenergetic (at indicated values of the Jacobi constant Cj) initial 
conditions taken on the x — y surface of section {px = 0, y > 0) for the planar circular (a— c), elliptic Sun- Jupiter (d— f) and and highly 
elliptic extrasolar (g— i) systems. In the elliptic case, initial conditions were computed with initial / = 7r/2. FLI were measured at the 
cut-off time of Tcut = 200 years. Very short-lived (scattering) trajectories that did not survive inside the Hill sphere for the cut-off time 
were discarded and are not shown. 



chaotic orbits with short (200 years in this example) resi- 
dence times inside the Hill sphere. These short-lived trajec- 
tories are the main contributions to the rapidly increasing 
number of chaotic orbits observed in Fig. 3] for e > 0.4. 
This does not mean, however, that capture will necessarily 
be enhanced, because a high density of strongly chaotic or- 
bits does not correlate with the number of very long-lived 
trajectories trapped close to KAM tori. 

On the other hand, Lyapunov exponents computed for 
individual trajectories cannot serve as a predictor of global 
stability, i.e. as an indicator of whether an orbit will stay 
long enough in a bound region, or if it escapes quickly. 
This Eispect in computing Lyapunov exponents for escap- 
ing (captured) trajectories concerns the notion of 'stable 



chaos' llMilani fc Nobilll992^ . This term was coined to refer 
to chaotic motions that are locally hyperbolic, but demon- 
strate macroscopic long-term stability, so that the Lyapunov 
time Tl (the inverse of the largest Lyapunov exponent Q) 
is substantially less than the 'event' time Te, e.g. the time 
for a trajectory to leave a certain region. However, there are 
indications of a possibl e simple linear correlation betwee n 
Te and Tl as shown in iLecar. Frankhn fc MurisonI Jl992l) . 
which suggests that knowledge of the Lyapunov time could 
allow one to predict an 'event' time. Concerning the resi- 
dence (escape) time near the planet, simulations in the pla- 
nar CRTBP, however, show that such a correlation with the 
Lyapunov time does not exist (Fig. El see also discussion by 
iMorbidelli fc Froeschla .1996i and .Varvoelis fc Anastasiadia 
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Il996f) . Among the randomly chosen escaping chaotic tra- 
jectories originating in the Hill sphere, there are many very 
long-lived examples that have diverse, even quite small, Lya- 
punov times. This means that local (microscopic) instabil- 
ities in the chaotic layer alone cannot explain the distribu- 
tions of survival probability. Rat her, as was pointed out by 
iTsiganis. Anastasiadis fc Varvog lis (2000), long-term trap- 
ping with short Tl near 'sticky' KAM structures can be at- 
tributed to the existence of phase space quasi-barriers and 
approximate integrals of motion. In the case of RTBP, these 
are the phase space structures around the saddle points 
at Li and L2 and the z-component of angular momentum 
with respect to the planet , respectively. In 3 D, the latter is 
approximately conserved dContopolouslllQGsl) for long-lived 
chaotic orbits which explains the unexpected strong corre- 
lation between initial and final inclinations of captured par- 
ticles ('inclination memorv'. .Astakhov et al...2003.'l . 



4.2 Monte Carlo simulations in the spatial 
ERTBP 

We simulated capture statistics in the spatial (three- 
dimensional) ERTBP by integrating isotropic fluxes of test 
particles that bombard the capture zone producing equal 
probabilities of initial conditions. Initial conditions were gen- 
erated as follows: the particle's position vector was chosen 
uniformly and randomly on the surface of the Hill sphere. 
Velocities were also chosen uniformly and randomly in accor- 
dance with the value of the CRTBP energy Ec- In turn, the 
starting energy (Jacobi constant Cj = —2Ec) was chosen 
uniformly random between its minimum possible value (as 
defined by the energy of Li when / = 7r/2) and its highest 
value (determined empirically such that above it the capture 
probability was essentially zero.) Then, by randomizing the 
initial true anomaly / we model equal chances for a test par- 
ticle to have any phase with respect to the mutual revolution 
of the primaries. The trajectories were integrated until one 
of the following occurred: the particle exited the Hill sphere; 
it penetrated a sphere, centred on the planet, of a given ra- 
dius (see Fig. |7] caption) ; it survived for a predetermined 
cut-off time. 

Choosing particles on the Hill sphere, where the mo- 
tion is chaotic or scattering, minimizes (although does not 
eliminate) the risk of acci dentally starting orbits inside im- 
penetrable KAM regions iNeto et alj|2004) . Although per- 
manently bound, these orbits could never actually have been 
captured because KAM regions cannot be penetrated at all 
in 2D and only exponentially slowly in 3D. It is only in the 
chaotic layer between scattering and stability that capture 
can happen. An initial swarm of incoming particles produces 
broad distribution of survivors with different residence times 
inside the Hill sphere. We monitored its dynamics up to cer- 
tain cut-off times (of the order of several thousand years) to 
find those long-lived chaotic orbits that may have been vul- 
nerable to capture by a relatively weak dissipative force such 
as gas drag. Only these orbits, in the CAC model, could be 
the precursors to the currently observed distributions of ir- 
regular satellites. In a statistical sense, any short-time fiybys 
are unlikely to have contributed to the primordial popula- 
tion of potential moons. Also, since the Sun- Jupiter system 
was chosen as an example, we eliminated test particles that 
penetrated the inner region of the Hill sphere occupied by 
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Figure 6. Escape time versus Lyapunov time for 2D CRTBP 
chaotic trajectories with initial conditions chosen randomly inside 
the Hill sphere. Escape was defined as the first crossing of the Hill 
sphere. 



Jupiter's most infiuencial regular moon Callisto. The mas- 
sive regular moons may have acted as a source of strong 
perturbation removing some temporarily captured moonlets 
following prograde orbits. This provides a possible reason for 
the observed prograde -poor distribution of jovian irregulars 
J Astakhov et alJl2003^ . 

Fig. shows the results of these simulations. The cap- 
ture probability density is plotted on the plane of initial 
and final inclinations. Unlike other properties, e.g., energy, 
inclination can be defined consistently in both the CRTBP 
and the ERTBP (e = and e > in isotropically pulsating 
coordinates). As the distributions in Fig|7] confirm, inclina- 
tion is approximately conserved during temporary capture 
in the ERTBP, so that the RTBP energy /inclination de- 
pendent CAC mechanism survives additional perturbations 
caused by deviations of the planet's orbit from circularity. 
This is clearly seen upon comparing Fig|71(a) and (b). The 
latter shows the capture probability for the Sun- Jupiter sys- 
tem with its actual eccentricity being used. We also ver- 
ified by similar Monte Carlo runs that the basic energy- 
inclination correlation discussed in Sec. 3 remains valid in 
the 3D ERTBP. We note, however, that since the energy is 
not a conserved quantity in the elliptic problem, a better 
representation of results is given in terms of inclinations. 

The most prominent effect observed in the ERTPB com- 
pared to the circular case is the overall decrease of cap- 
ture probability as eccentricity increases. Simulations indi- 
cate that, given the same time scale, the actual ellipticity 
of Jupiter's orbit accounts for approximately an order of 
magnitude lower capture probability than predicted by the 
circular model, while the relative number of progrades ver- 
sus retrogrades remains unaffected. 

Capture by typical extrasolar planets will be expected 
to become even more suppressed due to the wildly eccentric 
nature of these, essentially ERTBP, star-planet systems. To 
test this we used an eccentricity e — 0.4839 which is ten 
times greater than that of Jupiter's orbit but quite similar 
to what has been estimated, for instance, for several already 
detected extrasolars (the values of e are given in parenthe- 
ses): HD 108147 (0.498), HD 168443b (0.53), HD 82943c 
(0.54), HD 142415 (0.5), HD 4203 (0.46), HD 210277 (0.46), 
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HD 147513 (0.52), HD 190228 (0.5), HD 50554 (0.5), HD 
33636 (0.53). For consistency, we left unchanged all ol the 
other parameters and environmental factors (including Cal- 
listo whose hypothetical extrasolar analogs may well exist 
around exo-Jupiters) exactly as they were used in the sim- 
ulations for Sun- Jupiter. Fig|7Jc) indicates that, due to the 
ellipticity of the orbits and assuming similar time scales and 
other conditions, capture by extrasolar planets may be cis 
much as aproximately ten times less efficient than it could 
be by the giant planets of our system. 

The relatively low capture probability, predicted in 
ERTBP, adds to other destructive mechanisms possibly op- 
erating on as yet undiscovered irregular satellites of ex- 
trasolar planets. Th e loss of satellites through Yarkov sky 
l|Burns fc Cukll2002h or tidal (|Barnes fc O'Brienll200 ef- 
fects may also diminish the possibility for highly eUiptic ex- 
trasolar captors to marshal large populations of irregular 
moons. On the other hand, exo-planets with masses from 
one-half to four times Jupiter's mass on low eccentricity or- 
bits (e < 0.02) may be considered candidates for having 
irregular satellites. These may include, e.g., already known 
extrasolars such as HD 179949, 55 Cnc b, HD 169830c, HD 
187123, Tau Boo, HD 75289, 51 Peg, Ups And b, HD 195019. 



5 CONCLUSIONS 

The complexity of many-body systems evidently goes far 
beyond that of the simplest non-trivial case, the three-body 
problem. However, when considered on a hierarchical level, 
interactions between gravitional centers can often be de- 
composed to relatively low degree-of-freedom subsystems 
for which analysis from a dynamical systems point of view 
becomes possible. In the problem described here, the exis- 
tence of 'sticky' volumes of phase space where regular re- 
gions are surrounded by chaotic layers explains why several 
bodies may become temporarily (but for rather long time 
periods) trapped by mutual chaos-assisted capture. These 
quasi-stable configurations can subsequently be stabilised 
(or destroyed) on much longer time scales by various forces 
(dissipative or not). Here we have considered the effect in- 
troduced by a slow (compared to time-scales of motion in 
the chaotic layer) periodic parametric time dependence in 
the three-body problem. Our simulations show that chaos- 
assisted capture applies and is, in fact, only slightly per- 
turbed in cases when the primaries move on an elliptic orbit. 
It is expected to be of significant interest (e.g. in applications 
to essentially many-particle systems such as star clusters and 
the asteroid belt) to consider chaos-assisted capture at the 
next hierarchical level, i.e. when the primary two-body con- 
figuration is not restricted to a regular orbit, which may 
be the case in three-body encounters of comparably sized 
objects. 
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